Why a local moment induces an antiferromagnetic ordering: An RVB picture 
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Based on a Gutzwiller projected BCS wavefunction, it is shown that a local 5 = 1/2 moment is 
present around a vacancy site (zinc impurity) in a form of staggered magnetic moments, which is 
a direct consequence of the short-ranged resonating-valence-bond (RVB) pairing in the spin back- 
ground. 

PACS numbers: 74.20.Mn,74.25.Ha,75.20.Hr 



Experimentally, it has been observed that a zinc im- 
purity in the high- Tc cuprates generally induces a free 
magnetic moment in its neighborhood, which contributes 
to an additional Curie- like 1/T signature in NMR and 
NQR data In particular, the distribution 

of such a free moment around the impurity is in the form 
of staggered (AF) magnetic moments at the copper sites. 
Namely, if the free moment is polarized, say, by an ex- 
ternal magnetic field, its spatial distribution surrounding 
the zinc atom looks like an antiferromagnetic ordering of 
spins, with the total sum equal to the magnitude of the 
free moment. 

The induction of a local moment by a nonmagnetic 
zinc impurity illustrates the strongly correlated nature 
of the pure system and poses a great challenge to various 
high-Tc theories. Here two issues are involved. One is 
why a free moment should be spontaneously induced by 
a nonmagnetic zinc impurity. The second is that given 
the fact that a free moment is trapped nearby, why its 
distribution exhibits a staggered profile around the zinc 
impurity. 

There have been a number of articles on the micro- 
scopic origin of the zinc-induced local moment in the 
framework of the doped Mott insulator approaches 0, 0, 
1^. Wang and Lee [3 and Liang and Lee 8] attribute 
the induced local moment to a precursor effect of the AF 
ordering instability in the pure system. They consider 
the antiferromagnetic or spin density wave fluctuation in 
an RVB state nearby the zinc impurity. In this picture, 
the local moment is described by a local spin density 
wave order parameter which provides an explanation of 
the staggered profile of the local moment 13| . Qi and 
Weng p used a bosonic RVB theory to study a sud- 
den approximation for the zinc impurity state, in which 
an electron is annihilated from the zinc-free RVB state. 
In their picture, the staggered distribution of the local 
moment induced by the zinc impurity is a result of the 
short range AF correlations embedded in the pure sys- 
tem. They also argued that the phase string effect may 



protect the local moment for a topological reason. 

In this paper, we model the Zn-impurity by a 
Gutzwiller projected d-wave BCS wave function |0, ^ 
in the sudden approximation, in which the zinc doped 
system is constructed simply by locally removing an elec- 
tron at the zinc site from the pure system. We show that 
a free moment trapped nearby a Zn- impurity will gener- 
ally exhibit a staggered AF distribution, whereas the un- 
projected wave function does not. We show an explicit re- 
lation between the magnetic moment in the Zn-impurity 
state and the spin-spin correlation in the impurity-free 
RVB state. We have also examined a more sophisticated 
Zn-impurity RVB state, in which the BCS mean field 
state before the Gutzwiller projection is modified with 
some relaxation of the surrounding RVB background in 
the presence of an impurity (vacancy) site. The result is 
in further support of the sudden approximation. There- 
fore, the local staggered moment in the Zn-impurity state 
may be viewed as a consequence of the staggered spatial 
spin-spin correlation functions in the Zn-impurity free 
state. From this point of view, the staggered magnetic 
moments observed in experiments may not be necessarily 
related to the AF instability. Since the Gutzwiller pro- 
jection is essential to the staggered spin-spin correlation 
function, the Zn-impurity effect is a direct probe of the 
intrinsic RVB state in the cuprate superconductors. 

We start by approximately representing the pure sys- 
tem by a Gutzwiller projected d-wave BCS state 
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where Pq — J^j(l — nii^riii) enforces the no double oc- 
cupancy constraint with riia- denoting the electron num- 
ber at site z in a two-dimensional square lattice. The 
Gutzwiller projection transforms the singlet Cooper pairs 
in |BCS) into neutral spin RVB pairs at half-filling which 
persist into finite hole doping. Here the mean field state 
BCS) is derived from a renormalized mean field Hamil- 
tonian such as the Hamiltonian given in Eq. 

(4) below. 
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Now we construct a new state |5'o)zn removing lo- 
cally an electron on site io from |\l/o)- This state may 
be regarded as a "sudden approximation" for the true 
ground state |^)2n when site «o is occupied by a zinc 
impurity, similar to a construction previously used in 
the bosonic RVB approachj^. We note that |^o)zn 
the right spin quantum number (5* = 1/2) to describe 
a local moment. At the same time, both |vE')2n a-^id 
l*o)zn should reduce to |^o) in a distance sufficiently 
far away from the zinc site, as a single impurity is not 
expected to change the thermodynamic bulk properties 
of the system. Thus, we expect l'^)^^ and |5'o)zn should 
have a finite overlap such that the local states of them 
around the zinc impurity can be adiabatically connected 
by some kind of local perturbations, provided the local 
moment is present in |^)2n- Namely, if the S = 1/2 mo- 
ment is indeed distributed around the zinc impurity in 
the true ground state, then the "sudden approximation" 
state |'I'o)zn should be a good approximation as the first 
step to describe the related physics. 

Without the loss of generality one may define 



(2) 



Since IxPo) is a spin singlet state, one has 5*^ = ^1/2 
in |^'o)zn, which forms a doublet with the state ct^i |4'o)- 
Now let us inspect the distribution of this moment on the 
lattice by considering the quantity {S^)zn, the average 
of 5*1 in the state |4'o)zn- (Q) the average of 
operator Q in the pure RVB state jvPo), we have for the 
translational and rotational invariant RVB state. 



description of the t-J model, where the parameters t' 
and Aij can be determined by a set of self-consistent 
equations[ll[ll[l3. The results are shown in Fig. 1, 
which illustrates two things. First, there does exist stag- 
gered magnetic moments around the zinc impurity, which 
are actually determined by the short-range AF corre- 
lations already presented in the pure system as shown 
by the right-hand side of Eq.||3Jl; Second, the total mo- 
ment as the sum of these staggered moments is equal to 
— —1/2, according to the above construction. Namely 
the staggered moments spatially distributed around a 
zinc impurity and a "free" S* = 1/2 total moment can 
be naturally reconciled within this scheme. 
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FIG. 1: The spin distribution of {Sj'j around the zinc site 
(denoted by the plaquette in the center) in the state |'I'o)2n- 
The calculation is done on a 12 x 12 lattice of 22 holes, with 



= 0.1 and fr 



-0.2. 
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where n is the average number of electron per site. 
Therefore, the moment distribution is determined by the 
equal-time spin-spin correlation in the pure system. The 
right-hand side of Eq. (|3Jl can be numerically determined 
by a variational Monte Carlo calculation based on the 
Gutzwiller projected BCS state. In a system with anti- 
ferromagnetic spin-spin coupling on a square lattice, the 
sign of the spin-spin correlation is alternating in space. 
Therefore we expect a staggered magnetic moment in the 
Zn-impurity state. 

To make our discussion more quantitatively, we con- 
sider a renormalized Hamiltonian, from which the mean 
field ground state of |BCS) in Eq. (1) is derived, 



H 



MF 



(4) 



in which Ay denotes the d-wave pairing order parame- 
ter. This Hamiltonian can be taken as the mean field 



Physically, the local moment around the zinc impurity 
can be understood as originated from a broken RVB pair, 
with its one spinon partner (sitting at the site io) being 
removed as the consequence of the Zn substitution and 
the unpaired partner left behind to contribute to such 
a free local moment. Here the Mott physics also plays 
a crucial role for making such a free moment distribut- 
ing in the form of enhanced staggered moments as shown 
in Fig. 1. By contrast, without the Gutzwiller projec- 
tion operator Pq, one easily finds that the distribution 
of magnetic moments is given by 

(BCS|c|^Jg;)c,,T|BCS) 
(BCS|ct^^c,„T|BCS) 

|(BCS|cl^^c]jBCS)|2 + |(BCS|c1„^c,t|BCS)P 



< 



(5) 



which always remains negative and does not show a 
strong staggered oscillation like in Fig. 1. So the pres- 
ence of staggered AF moments around the zinc impurity 
is a very important feature for an RVB superconductor 
(Gutzwiller projected BCS state) |5'o), but not for a con- 
ventional BCS state |BCS). 
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The above discussion is based on the assumption that 
l^o)zn ^ finite overlap with the true ground state 
4')2jj- Namely, the existence of a local moment is as- 
sumed to be around a zinc impurity. If the local moment 
in |'I'o)zn' created by the "sudden approximation" from 
the original ground state |^'o), eventually escapes to in- 
finity (boundary) in the ground state, then |^'o)2,i and 
h^')zn ^^'c orthogonal and the present "sudden approxi- 
mation" used above will fail. In principle, an RVB state 
does not prevent such a local moment (spinon) from es- 
caping away to the boundary as it can re-arrange the 
RVB pairings in the bulk. So a more sophisticated ap- 
proach is needed in order to further demonstrate the sta- 
bility of the local moment around the Zn impurity. 

For this purpose, in the following we consider a mod- 
ified wave function for describing the local moment. In- 
stead of removing an electron in the projected state of 
the pure system, the impurity state is constructed as 
a Gutzwiller projected state in which the unprojected 
mean field state is obtained in the presence of the zinc 
impurity with a local quasiparticle excitation. Generally, 
an s-wave unitary scatter doped into a d-wave supercon- 
ductor will induce a low-energy d-wave resonance state 
around it, ensured by the d-wave pairing svmmetrv|l5|. 
and if we denote the creation operator of this resonant 
state as 7]., then the modified wave function for describ- 
ing the local moment state is given by Pelt \^^^)ziv 
which |BCS)2n denotes the mean field ground state in 
the presence of the zinc impurity. Note that here the 
occupied d-wave resonance state will have the right spin 
quantum number to describe a local moment and at the 
same time any double occupancy of the same local state, 
which would destroy the moment, is prohibited by the 
Gutzwiller projection. The remaining task is to check if 
the proposed wave function still exhibits a staggered dis- 
tribution of the moment. The mean field description of 
the single zinc impurity at site io is given by Hmf of Eq. 
(4) with the impurity site iq excluded. 
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FIG. 2: The distribution of (Sj) around the zinc site in the 

state PgtI |BCS)2„. The parameters used are the same as in 
Fig. 1. 



The mean field ground state of H^p and the d-wave 
resonant state 7^. can be obtained numerically on the fi- 
nite lattice. Then one can perform the VMC calculation 
to implement the Gutzwiller projection which is done on 
a 12 X 12 lattice with 22 holes. Fig. 2 shows the spin 
distribution around the zinc site in Pg7| |BCS)2„- The 
spin distribution in Pg7| |BCS)2jj is slightly more ex- 
tended than that in |^o}z,i as a result of the relaxation 
of the RVB background but otherwise remains quite sim- 
ilar. Since both Pg7| |BCS)2„ and |5'o}zn have the same 
spin quantum number and show similar spin distribution, 
it is natural to expect that they are not very different 
from each other. In fact, we have calculated their over- 
lap by using the VMC and the result is shown in Fig. 3. 
From the figure one sees that the overlap actually satu- 
rates to a finite value about 0.8 in the thermodynamic 
limit. Thus the two states are indeed very close to each 
other in terms of both the large wave function overlap as 
well as the staggered profile of the local moment. 
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FIG. 3; The overlap between Pg7| IBCS)2„ and |*o)zn as a 
function of the inverse lattice size . The calculation is 
done with the hole density kept approximately at 15 percent 
(for example, 62 holes in the 20 x 20 lattice). The parameters, 
jj- and jr, are the same as in Fig. 1. 

Energetically a close competitor of the state 
Pg7| |BCS)2„ is the projected mean field ground state 
Pg |BCS)2„. Since the mean field excitation energy of the 
resonant state created by 7j is small, we expect that the 

energy difference between Pg7| |BCS)2n and Pg |BCS)2„ 
be small also. Indeed the energy difference in the VMC 
calculation based on the t — J model is generally very 
small, although the energy of Pg |BCS)2„ remains to be 
lower. So strictly speaking, how the local moment can 
be trapped near the zinc impurity, due to the distortion 
of the RVB order parameter around the impurity is still 
unsolved in the present approach. In contrast, a topolog- 
ical reason coming from the charge degrees of freedom in 
the bosonic RVB theory provides the protection of the 
local moment by responding nonlocally to the presence 
of a vacancy in the bulk. Whether the underlying physics 
for the trapping of the local moment is due to the local 
RVB distortion or some more subtle topological reason 
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needs to be further studied in the present fermionic RVB 
framework. Nevertheless, once a local moment is trapped 
around the zinc impurity, a staggered moment profile al- 
ways emerges, as a direct consequence of the presence of 
RVB correlations in the pure system. 

In summary, we studied two fermionic RVB wavefunc- 
tions describing the zinc doping-induced local moment 
effect in the doped Mott insulator. We find that the 
AF staggered profile of the local moment distribution 
emerges naturally and can be generally attributed to the 
short-range AF correlations of the RVB structure in the 
pure system. According to this description, the existence 
of the induced local moment with a staggered profile in 
the cuprates provides a strong experimental indication 
for the underlying RVB structure in the zinc-free sys- 
tems. 
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